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Quantum decoherence in the entanglement entropy of a composite particle and its 
relationship to coarse graining in the Husimi function 
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I investigate quantum decoherence in a one-body density matrix of a composite particle consisting 
of two correlated particles. Because of a two-body correlation in the composite particle, quantum 
decoherence occurs in the one-body density matrix that has been reduced from the two-body density 
matrix. As the delocalization of the distribution of the composite particle grows, the entanglement 
entropy increases, and the system can be well-described by a semi-classical approximation, wherein 
the center position of the composite particle can be regarded as a classical coordinate. I connect the 
quantum decoherence in the one-body density matrix of a composite particle to the coarse graining 
in a phase space distribution function of a single particle and associate it with the Husimi function. 


In recent decades, quantum entanglement has at¬ 
tracted a great deal of interest in various fields. To es¬ 
timate correlations in quantum systems, entanglement 
measures such as the entanglement entropy (EE) have 
been intensively studied lj-flj. Many entanglement mea¬ 
sures are defined by reduced density matrices which de¬ 
scribe the structure of the Schmidt decomposition and 
contain information about entanglements in quantum 
systems. In entangled states, the EE is produced by 
quantum decoherence caused by a reduction in the num¬ 
ber of degrees of freedom (DOF). In my previous papers, 
I calculated entanglement measures of the one-body den¬ 
sity matrix in nuclear systems mini and showed that 
the EE is enhanced by the delocalization of the distribu¬ 
tion of clusters, which are composite particles of spatially 
correlated nucleons. My first aim in the present paper 
is to understand how quantum decoherence occurs and 
how the EE is produced in the one-body density matrix 
of correlating particles. 

Quantum decoherence—that is, the quantum 
entropy—has been also investigated with the coarse 
graining of distribution functions in a phase space. 
The Husimi function fl2l4l4l ] is known to have finite 
the Wehrl and Renyi-Wehrl entropies [HI, HU because 
of the coarse graining by a Gaussian smearing of the 
Wigner function. It has been shown that the Wehrl and 
Renyi-Wehrl entropies are increased by delocalization of 
the distributions in quantum systems [l7l - l20| . Campos 
et al. have discussed a correlation between the EE and 
the Renyi-Wehrl entropy in entangled states j2c}. One 
of the fundamental questions in quantum physics is how 
the quantum decoherence in the reduced density matrix 
of entangled states can be connected to the coarse 
graining of distribution functions. My second aim in this 
paper is to understand the correspondence between the 
quantum decoherence in the one-body density matrix 
of correlating particles and the coarse graining in the 
Husimi function of a single-particle state. 

In this paper, I investigate the EE of the one-body den¬ 
sity matrix of a two-body system in which two particles 
are strongly correlated to form a composite particle, and 
I discuss how quantum decohence occurs in the reduc¬ 
tion of the DOF. To describe two-body wave functions, 


I adopt a cluster wave function in the generator coordi¬ 
nate method in nuclear physics [2l], H3 ■ Let us consider 
a system where two particles (cl and c2) with masses 
to and um form a bound state with an attractive inter¬ 
particle force. I assume that the bound state is described 
by the lowest state of a harmonic oscillator (ho) potential 
and can be approximately treated as an inert composite 
particle, where intrinsic excitations cost a relatively high 
amount of energy compared with the center of mass (cm) 
motion of the composite particle. In this approximation, 
a total two-body wave function is given as 

\^) = f dsF( S )\s-,b) 1 \s-,b 2 ) 2 , (1) 

(ri\s;b)i = /(b 2 tt) 1/4 , (2) 

where & 2 = b/y/u , and (VRl 2 ) |\I>( 2 )) = 1. ri(r 2 ) is the co¬ 
ordinate of cl(c2). Here, I describe the one-dimensional 
case, but the present model can also be extended to the 
three-dimensional case. |s; 6)i|s; 62)2 indicates the com¬ 
posite particle localized around the mean position s, and 
\j/(2) gj ven by the superposition of different s states 
with the weight factor F{s). I^ 2 )) can be expressed by 
the cm motion and the intrinsic wave functions as 

l* (2) > = |<Mfl)>|& nt (r)>, (3) 

(R\*g(R)) = J ( 4 ) 

with the cm coordinate R = u\r\ + ti 2 r 2 {u\ = l/(w + 1) 
and U 2 = u/(u + 1)), the relative coordinate r = rq — r 2 , 

and b G = y/ufb. Here, (r|(/)i nt (r)) = e ^ /(bln) 1 ^ 
with b r = b/yfuJ, is the lowest intrinsic state for the ho 
potential, Uh_o{ti,b r ',r) = —h 2 r 2 /2[j,bf,, with [i = u 2 ?u. 
Thus, general low momentum states of the inert compos¬ 
ite particle can be expressed by the form a, in which the 
cm motion 4 >g( 1?) is expressed by the shifted Gaussian 
expansion as given in Eq. a. 

The one-body density matrix p^ (2) for cl is defined by 
the matrix reduced from the many-body density matrix 
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P%) = |* ( 2 ) )(tf (2) | as p$ 2) = Tr 2 [p^ 2) 2) ] and is given as 


P%) = J dsds'W(s',s)(s';b 2 \s- 1 b 2 )\s;b}(s';b\, (5) 
= (<h\P%) Wi) 


(i& 2 7r ) 1 / 2 




( 6 ) 


where W(s',s) = F*(s')F(s) and Trp^ (2) = 1. Note 

that p^( 2 ) for u = 3 equals to the one-body density ma¬ 
trix of an a cluster composed of four nucleons with an 
equal mass investigated in previous papers US El- The 
Wigner transformation (Wigner function) of p^\ 2 ) is 

P W (p*m]Qi,Pi) = ( dr)(qi + ^\p^L\Qi - 


J q/ ( 2 ) i 


xe 


-^■(9 1 - s ) 2 -5p-(9 1 -s') 2 --p-{p i -^-(s-s')} 2 


(7) 


The Renyi EE of order 2 (Renyi-2 EE) and von Neumann 
EE for 'fh 2 ) with the one-body density matrix p^\ 2) are 
given as 


S™(p%) = -\n. (Tr^l,} 2 ]) 

= -In i^J dqidp 1 {p W (p ( ^ 2) ;q 1 ,p 1 )} : 


S VN (P%) = -Tr 


o (1) lno (1) 


( 8 ) 

(9) 


If p W (/3^( 2) ; qi,Pi) > 0 is satisfied in the entire phase 
space, I can consider the phase-space Shannon entropy 
S sh (p{q,p)) = f ^£p(q,p)lnp(q,p) for p w (p$ 2) ; gi,Pi) 
as, 


F(s') ss F(s). Then I obtain a semi-classical approxima¬ 
tion, 



where f(s) oc F(s) whose normalization is determined by 
/ dq 1 p^ 2 f(q 1 ,qi) = / ds|/(s )| 2 = 1. This corresponds 
to 


P%) « P l Sif = J ds\f(s)\ 2 \s;b r )(s;b r \, ( 12 ) 


where the parameter s and |/(s )| 2 are regarded as a clas¬ 
sical coordinate and a classical distribution of the com¬ 
posite particle, respectively. In the large u limit—that 
is, the large c2 mass limit— p(qi,q[) P cl (qi,q'i ) and 
b r —»• b. Note that, |s) and |s'} are not orthogonal to 
each other because of the quantum fluctuations of the cl 
position in the composite particle. 

In the semi-classical approximation given by Eq. m, 
the Wigner function is approximated as 


P {p^m,qi,Pi) 


ds\f(s)\ 2 e 


1 / \2 b r 2 

— 7t( s_ 9i) ~tzPi 


, (13) 


which is nonnegative definite. Using p^ (2) cl I also define 
the EEs in the semi-classical approximation: S ' 112 ’ 01 = 

S R2 (P ( X) and ^ W - Sh > cl = 5 W - Sh (p^f) ; 

For simple examples, I first consider the zero- 
momentum state of the composite particle in a finite 
volume V described by a constant F(s). I assume that 
b >6 and the contribution of the box boundary can be 
ignored and obtain 


^( 2 ) (91,^) = 
P KPismWuPi) = 


1 -l^0n-9i) 2 

V 

2&J.7T 1 / 2 2 

— y —e ^ Pl 


(14) 

(15) 


5 W - Sh (p%) = S Sh (p w (p^( 2) ; ?i,Pi)^ , (10) 

which I call the “Wigner-Shannon EE”. 

In the one-body density matrix p^ (2) and its Wigner 
transformation, quantum decoherence occurs and pro¬ 
duces the EEs because of the factor (s'■, b 2 \s- : b 2 ) = 
exp[— j§?(s — s') 2 ], which originates in the reduction of 
the DOF of c2. Indeed, in the case of u = 0, without this 
factor, pT) = {pT )} 2 and the Renyi-2 and von Neumann 
EEs are zero that corresponds to a pure single-particle 
state. 

Let us consider a semi-classical approximation of p^ (2) . 
The factor exp [—^(s — s') 2 ], which is the source of 
the quantum decoherence, has a sharp peak around 
s' ~ s with a width 2 b/y/u. I assume that the 
function F(s ) is a slowly varying function compared 
with exp [— 45 ?(s — s') 2 ], arL( l it can be approximated as 


In this case, p^} 3) (qi , q [) = p^}if{qi, q[) is satisfied. The 
Renyi-2, von Neumann, and Wigner-Shannon EEs are 

S R2 = lnU eff - iln(27r), (16) 

S vN = 5 w-Sh = S R 2 + l( 1 _ ln2 ) ! (1.7) 

where Kff = V/b r denotes the effective Volume size for 
the cm motion. These results are not valid for a small 
V eS because of the box boundary. 

The one-body density matrix is diagonalized in the mo- 

mentum space with a Gaussian distribution, exp[— jfip 2 \- 
This indicates that the one-body density matrix of a free 
composite particle is equivalent to a thermal state of a 
single particle at finite temperature kT = H 2 /2mb 2 . The 
temperature is of the same order as the mean kinetic en¬ 
ergy, h 2 /2mb 2 1 of constituent particles confined in the 
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composite particle. Strictly speaking EEs are not ther¬ 
modynamic entropies, however, by associating the one- 
body density matrix of the free composite particle with 
a quantum mixed state of a single particle, I can propose 
an interpretation of the entropy production and thermal- 
ization as follows: when the DOF of c2 are reduced, the 
quantum decoherence occurs, producing the entropy, and 
simultaneously, the intrinsic kinetic energy of the com¬ 
posite particle is converted into heat. 




FIG. 1: EE of a composite particle with a Gaussian distribu¬ 
tion for (a) u = 1 and (b) u = 8 . e s for the Renyi-2 (S R2 and 
S' R 2 ’Ci), von Neumann (S vN ), and Wigner-Shannon gW -311 ) 
EEs are plotted as functions of the effective volume size, v e ff- 
The Renyi-2 EE in the large u limit is also shown in the panel 
(b). 


Next, I consider a composite particle moving in an ex¬ 
ternal ho potential, where the lowest state of the com¬ 
posite particle is given by the Gaussian distribution, 

F(s) = e~ zb 1 / ( B 2 7 t) 1 / 4 . This gives the exact solution to 
the two-body wave function, (n, = (i?, rjH/G)), 

for ho potentials Uh 0 (M, /?; R ) + Ug 0 {p, b r ; r), with M = 
( u + 1 )m and /3 = \/ B 2 + u\b 2 . In the B = 0 limit, 
describes a localized composite particle that corre¬ 
sponds to a non-entangled (uncorrelated) state of two 
constituent particles and has zero Renyi-2 and von Neu¬ 
mann EEs. As B enlarges and the delocalization of the 
cm of the composite particle grows, the EEs increase. 
The Wigner function and EEs for are 


P W (/S*(2);9i’Pi) 

S R2 

gW-Sh 


VWTb 2 

iln(l + w 2 ff ) - ilny, 
S R2 + 1 - ln2, 


(18) 

(19) 

( 20 ) 


where 7 = (1 + (u + l)u 2 ff )/(l + uv 2 s ), and v e g = B/b 
denotes the effective volume size. The EEs increase as 
v e g enlarges and approaches lnu e ff in the large v e g limit. 
In the semi-classical approximation, the Wigner function 


and EEs are 


9 W (b (1) ■ 


<?1, Pi) 

cl 

2 b r 


sjbl + B 2 

S R2 

cl 

oR2,cl _ 1 
2 

gW-Sh 

cl 

gW-Sh,cl _ 




, ( 21 ) 


where v c .eff = B/b r . I show the EEs for u = 1 and u = 8 
in Fig. [T] S' vN is calculated numerically, as was done 
in the previous paper 0 - S R2,c\ fol . the 

semi-classical 

approximation agrees well with S R2 in the > 2 case 
to within 10 % error for u = 1, and the agreement is 
better for the larger mass ratio, u = 8 . 5 ,w_Sh has a 
constant shift 1 — ln 2 (a constant scaling e /2 in the e s 
plot in Fig. |T]) from S^ 2 , and it is finite even at w e ff = 0. 
S vN starts from zero at w e ff = 0 and approaches S w ~ sh 
as v e ff increases. As the mass ratio u increases, the EEs 
converge on values in the large u limit. 

Finally, I connect the quantum decoherence in the one- 
body density matrix of the composite particle to coarse 
graining in the phase space distribution function of a sin¬ 
gle particle and associate it with the Husirni function. 
Let us start from the Wigner transformation of the full 
two-body density matrix (2) , 

P W (p%)',qi,Pi,Q2,P2) 

= J di^(qi + + |) 

xe -ipiv-iP2^ (24) 


It is rewritten by a separable form in the phase space for 
the cm and relative coordinates as 


P W iP < i(2)',qi,Pi ,q2,P2) 


p w V>%>;Q,p)p w (P% 


q,p), 

(25) 


where p ( #J = |4>g)('I > g| and p^ = |</> in t> | - The 

Wigner function of the one-body density matrix can be 
written as 


W/- 41 ) \ [ dq2dp2 w/-(2) \ 

p (Pym\q\,pi) = 97rh p {pym-,qi,Pi,q2,P2) 

~(r) Ql Q 


X P UVi 


U2 


,Pl ~ UiP 


-p^(9i-Q) 2 -§(pi-«i ^) 2 

xe r * 


(26) 


(27) 


Here I use relations q = (qi — Q)/u 2 , p = Pi — u\P, 
and the transformation dq^dp 2 = \J\dQdP with the de¬ 
terminant of Jacobian |J| = l/u 2 - This means that 
P W {iyl 2 ) \qi,Pi) * s regarded as a coarse grained distribu¬ 
tion function of p w (p$(?; Q 1 P) with a Gaussian smear¬ 
ing. In other words, the quantum decoherence caused 
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by the reduction of the DOF of c2 can be interpreted 
as the coarse graining in the phase space distribution of 
a single-particle state. It is important that, if the in¬ 
ternal DOF are decoupled from the cm motion of the 
composite particle, Eq. (l26l) describes a general form of 
the coarse grained distribution function that corresponds 
to the Wigner function of p^ {2) . 

Let us consider the u = 1 case and associate the coarse 
graining in Eq. (1271) with the Husirni function. Eq. (1271) 
for u = 1 is rewritten as 


P W {p%)\<h,Pi) = 2 Pft/ 2 (P^ 6 G;9i>JPi)> 
pr/2(p (1) ;«;9^) 


SWm!^¥'p- ¥■ 


iPr, 

HJ2 


xe 




(28) 


(29) 


where pf /2 is normalized as / ^f0TjPhl2 = L Pn /2 even - 
tually has the same form as the normal Husirni func¬ 
tion, except for scaling of Planck’s constant h —> h/2. 
I call “fi/2-Husirni function”. It is clear that 

P W (P^ 2 ) ;<7i,Pi) for the two-body state is equivalent 
to twice the ?i/2-Husimi function for the single-particle 
state, |'I>g). The Gaussian smearing in the coarse grain¬ 
ing originates from the reduction of the DOF of c2 in 
p w (pJ2 t ; ( q,p )) as shown previously. Note that the h/2- 
Husimi function is not a distribution function for a phys¬ 
ical single-particle state, but is regarded as a “distribu¬ 
tion” function defined in the down-scaled phase space, 
h —» H/2. The reason for the down scaling h —> h/2 is 
that the (< 72 , P 2 ) phase space is scaled down in the trans¬ 
formation from (qi,pi , 92 ,^ 2 ) to (qi,pi,Q,P). 

Considering the one-to-one correspondence between 
P W (PyU',QuPi) and P^ /2 (p < ^-,qi,Pi), I can connect en¬ 


tropies defined by the ft/2-Husimi function to EEs as 


yWehrl/ ~(G)\ 


> V 2 W G J ) = 5 w - bh (4%)) + ln2 , (30) 
S^ ehl \& = S R2 (p%) + ln2. (31) 

Here I dehne the Wehrl entropy and the Renyi-Wehrl 
entropy of the order 2 in the down-scaled phase space as 

5 We hrl (^ (1) ) = _ J Pn/2 ( P {1) l g. P) h Wn/2 (g» P) . 


^R2-Wehrl 
5 h/2 


(p (1) ) = -111 
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In summary, I investigate the quantum decoherence in 
the one-body density matrix of the composite particle 
that comprised two correlated particles in the inert com¬ 
posite particle approximation. Because of the two-body 
correlation in the composite particle, the quantum deco¬ 
herence occurs by the reduction of the DOF of the second 
particle. As the delocalization of the distribution of the 
composite particle grows, the entanglement entropy in¬ 
creases. I found a one-to-one correspondence between 
the quantum decoherence in the reduced density matrix 
and the coarse graining in the phase space distribution, 
which is related to the Husimi-like function defined in the 
down-scaled phase space. In the present paper, the in¬ 
ert composite particle approximation is applied to static 
systems but it can also be extended to time-dependent 
systems if the energy scale of the internal DOF of the 
composite particle is decoupled from that of the exter¬ 
nal DOF. The present study may shed light on the fun¬ 
damental problems of quantum decoherence and coarse 
graining which produces entropies in quantum systems. 
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